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TECHNICAL NOTE 277^ 


A METHOD FOR FINDING A LEAST-SQUARES POLYNOMIAL 
THAT PASSES THROUGH A SPECIFIED POINT 
WITH SPECIFIED DERIVATIVES 
By Neal Tetervin 

SUMMARY 


A method is presented for finding an mth.-degree polynomial that 
passes through a specified point with (m - j) specified derivatives 
(1 = J = m) and is a least-squares polynomial for other points spaced 
at unequal Intervals of the Independent variable. 


INTRODUCTION 


Recently it became necessary to be able to find the velocity and 
the rate of change of velocity along the surface at any point on the 
forward portion of a symmetrical airfoil at zero angle of attack. The 
given data consisted of values of the velocity at unequal intervals 
along the surface and of the first derivative of the velocity at the 
stagnation point. The given velocity at one of the points seemed to be 
subtly in error. 

In the development of a method for finding the velocity and its 
first derivative, a procedure was needed for obtaining a third-degree 
polynomial that passes through a specified point with a specified slope 
and is a least-squares curve for other points spaced at unequal intervals 
of the independent variable. A search of the available literature 
(refs. 1 to 3) in which the least-squares method is discussed failed, 
however, to disclose a procedure that ccanbines the least-squares method 
with the requirement that the polynomial pass throu^ a specified point 
with a specified first derivative. Such a procedure, therefore, had to 
be developed. 

In the course of the analysis, it became evident that its scope 
could easily be e 2 q>anded to provide a method for finding an mth-degree 
polynomial that passes throu^ a specified point with (m - j) specified 
derivatives (l ^ j = m) and, moreover, is a least-squares polynomial for 
the remaining points. The specified (m - j) derivatives are not neces- 
sarily the first (m - j) derivatives. 
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The method presented should he useful when the value of a function 
or the value of a function and some of its derivatives are known from 
theory at one value of the Independent variable and experimentally 
obtained values of the function are available for other values of the 
independent variable. An example is the idtting of a polynomial to the 
portion of a laminar-boundary-layer velocity profile near the surface. 

In this case the value of the velocity and some of the derivatives of 
the velocity are known at the surface from theoiy and measured values 
of the velocity are available at a number of points throu^ the boundary 
layer. The procedure should also be useful when it is necessary to 
divide the total range of the independent variable into convenient 
Intervals and to find a least- squares polynomial for the points in each 
interval. The polynomial can then be made to pass throu^ the point 
given by the preceding polynomial at the end of the preceding interval 
with the required derivatives. Examples are the fitting of an airfoil 
velocity distribution by a number of polynomials and the fitting of a 
boundary- layer velocity profile by polynomials. 


SYMBOLS 


N + 1 

number of points 

m 

degree of polynomial 

X 

independent variable 

y 

dependent variable 


specified point throu^ which curve passes 

Z = X - Xq 


1 

II 

a 


Co^Ci,C2, . . .C^ 

coefficients in polynomial 

^u^^^U2^^U2^ • * 

coefficients to be calculated by least-squares 
method 


coefficients to be calculated from given 
derivatives 

Cp 

pth coefficient in polyncmlal 
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m - j 

n = 0, 1, 2, 3 , . 
S 


number of coefficients to be calculated by 
least-squares method 

number of coefficients to be calculated from 
given derivatives 




sum 


ANALYSIS 


« 


Let there be given N + 1 points (xo,yo)j 

It is required to find an mth-degree polynomial that passes throu^ the 
point (xQ,yQ) with specified (m - j) derivatives and is also a least- 

squares polynomial for the remaining N points. Some of the 
remaining N points can lie on one side of Xq and some on the other 

side of Xq, or all the remaining N points can lie on one side of Xq. 
The polynomial is written in the fom 

y = Cq + CjZ + C2Z^ + + . . . + (l) 

where 

Z = X - Xq 


The requirement that the cuiTre determined by the polyncmial 
(eq. (l)) pass throu^ the point (xQ,yQ) makes Cq = jq . Additional 

(m - j) coefficients in the polynomial (eq. (l)) are found from the known 
(m - j) derivatives at point placing z = 0 in the relation 

for the pth derivative: 


^ = pICp + (p + 1)!C 3_z + 

dz-^ 


(p + 2^)’ O 2 (p + 3)! 
^+2 2 + 


21 


3! 




(p + O- 




m; 


(m - p)‘ 




m-p 


( 2 ) 
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Thus 



(3) 


There are thus (m - j + 1) known coefficients. The total number of coef- 
ficients in the polynomial (eq,. (1)) is m + 1; therefore^ j unknown 
coefficients remain. 

For the least-squares method to be applicable, the number of unknown 
coefficients, j, must be less than N. When j := N, the polynomial passes 
throu^ each point and has the specified derivatives at point (xQ,yQ^. In 

this case, the j unknown coefficients can be found more simply than by 
the least-squares method. In the present analysis, j < N. 

The j unknown coefficients are to have values that make the poly- 
nomial (eq. (1)) a least-squares polynomial for the N points 
(Xi,yi), unknown coefficients are 

to have the values that minimize the sum of the squares of the differences 
between the y values calculated by equation (1) and the given y values 
for Xj_, X 2 , . . . Xjjj thus. 


N / .2 

S = ZI(Co + C^Zn + C 2 z/ + + . . . + - yj (4) 

n=0 


is to be a minimum, where 

2n = ^ ^ 

The substitution of 0 for 1 can be made in the lower limit of the 
summation because Zq = 0 and Cq = yQ. The ejqpression for S can be 
written as 
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where C^ , , . . • C^ . are the j unknown coefficients and 

123 J 

Cq and Cjj. ^ are the (m - j + 1) known coef- 
ficients. The integers U 2 , used as subscripts and 

superscripts are not necessarily consecutive. Now let 


Equation (5) then becomes 


/ u-i Uo n-3 Ua \2 , 

S = _2 + CugZn. + Cu^Zn -^ + • . . + Cu^Zn ^ - w^j (7. 

The sum S is a function of the j unknown coefficients C^^, *^U2^ 

Cjj^ , . . . . For S to be a minimum it is necessary that 

3 J 

^ = 0, = 0, = 0, . . . -is_ := 0 (8] 

SCu2 ^^^^3 

The application of the conditions given by equations (8) to equation ( 7 ) 
results in the following equations: 

^ U-i Up Uo Uj \ U-| 

\p”l^^ ^ *^U2^ *^U2^n + • . . + C^^Zji ” ■‘^nj^n ” ® 

/ Ut Up Up U^ \ Up 

^(Cu]_zn + Cu2Zn + ^ + ■ ■ ■ + C^^Zn ^ - w^jz^ ^ = 0 

r 

^ ^^uiZn ^ + C'U2^n ^ + *^u^Zn 3 + . . . + Cu^z^ j - Wn^Zji 3 = 0 
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s i * =.2 i • s i: i ■ iz -a ^ 

■‘- n=0 n=0 ^ n.*0 J n=0 n=0 

Oui ± %“"'■"= * S i * “«3 i *••■-=., i " ± 

n=0 ^ n=0 n=0 J n=0 n=0 


+U^ — - V ■ Urt+U^ 2 u 

Cu, ^ Zn ^ 3 + Cuo H % ^ 3 + C„ y_ Zn 

i r^^-A ^ v»»^A J T\— A 


i ^"“3 . , . , * c„ t: C 3*“3 . i: ,„.„"3 

^ n=0 " n=0 n=0 


± • C i * S i - c^, i i ..O 


Now, for simplicity in writing, introduce the notation 

N e N 
ZI ^ = a, 
n =0 *' ^ 


21 ^n^n^ ^ ) 

n =0 ' 

Equations ( 10 ) foi ^u.2^ • ’ * ^Uj hecome 

^uj^'^2u]_ ■*■ ^U2®U2 +u^ ■'■ + . . . + ^Uj®Uj+Uj_ " ^U]j_ 


*^Ui'^Ui-fU2 ^U2°'2 u 2 *^U2'^U2+U2 + • • • + *^Uj'^Uj-t-U2 ~ 

*^Ui*^Ui-fa 3 ■'■ ^U 2 "^U 2 +U 3 *^U 3 '^ 2 u 3 + • • • + ~ ^u- 




^Ut °’l 


2 _°^U 2 _+Uj ^U 2 °^U 2 -t-Uj ■*■ + . . . +• Cu^^ 02 uj ~ 
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The values of C^l^, 



are determined, hy solving 


the J simultaneous linear equations (eqs. (ll)). A convenient method 
of solution is given in chapter I of reference 3« 


NUMERICAL EXAMPLE 


A simple example is given to illustrate the use of the method for 
a specific case. 

A thiitL-degree polynmial is to pass throu^ the point ^xq = 
y-Q = 125) with a slope of 75 and is to he a least-square polynomial 
for the points: 


X 

y 

1 

1 

2 

8 



7 

3^3 

8 

512 

9 

729 

10 

1000 


The polynomial is 


where 


y = Cq + Cq^z + C2Z^ + C^z 


The requirement that y = 125 


z = X - 5 

for X = 5 (z = 0) leads to 


Co = 125 

Fi;om equation (3) with p = 1 and ( — ) = 75^ the resulting value 

obtained for Cj^ is 


Cl = 75 
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In tills example there is one known derivative and the polynomial is of 
the third degree; therefore, 

m - J = 1 
m = 3 


Consequently, j = 2 and 



With Cq 


125 and 


Cq_ = 75 ^ equation ( 6 ) becomes 


^n = yn - 125 - 75 zn 


In order to sinpllfy the computation of C2 and C^, the following 
table was prepared: 


n 

% 

m 


z 3 
Zn 

zn^ 

Zn5 

z ^ 

Vn 

ynZn^ 

ynZn^ 

B 

5 

0 

0 

0 

0 

0 

0 

125 

0 

0 

B 

B 

D 

16 

-64 

256 

-1,024 

4,096 

m 

16 

-64 

2 

2 



-3 

9 

-27 

81 

-243 

729 

8 

72 

-216 

3 

B 

H 

B 

B 

1 

1 

H 

1 

64 

64 

-64 

B 

7 

2 

B 

8 

16 

32 

64 

343 

1,372 

2,744 

5 

8 

3 

9 

27 

81 

243 

729 

512 

4,608 

13,824 

6 

9 

4 

16 

64 

256 

1,024 

4,096 

729 

11, 664 

46 , 656 

7 

10 

5 

25 

125 

625 

3,125 

15,625 

1,000 

25,000 

125,000 

> 

1 

1 

80 

132 

1,316 

3,156 

2 ^, 3 kO 


42, 796 

187,880 
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With. u-]_ = 2, U2 = 3^ = 3 equations (ll) become 

Cgo^lj. + ^3^5 ~ ^2 
02^5 + C^cTg = 02 


where 



1,316 


= 21 ^rP = 3,156 

n=0 

<^5 = ^ = 25,3it0 

n=0 


and 






- 125 - T 5 zJz ^2 


7 _1 

y y z 2 _ 125 

•'n n 

n=0 


7 

75 y~ z 

n=0 


3 

n 


= 42,796 - 125 ( 80 ) - 75(132) 
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and 


i^3 = 51 ^n^n^ 

n=0 


= i (yn - 125 - 752n)%5 


n=0 


= f: yn^3 . 125 £; . t 5 i; 


n=0 


n=0 


n=0 


= 187,880 - 125(132) - 75(1316) 


= 72,680 


The equations for C2 and then are 


I316C2 + 3156C3 = 22,896 


3156C2 + 25,3^003 = 72,680 


Solving for C2 and C3 results in 


C2 = 15 


C3 = 1 


The required cubic equation then is 

y = 125 + 75z + 15z^ + 


or 


y = 125 + 75(x - 5} + 15(x - 5)^ + (x - 5)3 
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vhich reduces to 


y = 


x3 


the equation frcorn. which the set of values of y 


against x was obtained. 


CQHCLUSIONS 


A method is presented for finding an mth-degree polynomial that 
passes throu^ a specified point with (m - j) specified derivatives 
(1 ^ m) and is a least-squares polynomial for other points which are 

spaced at unequal intervals of the independent variable. 


Langley Aeronautical Laboratory 

National Advisoiy Committee for Aeronautics 
Langley Fields Va.^ June 1952 
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